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Abstract 

Using a numerical approach we examine the dynamics of zz correlations for the spin- 
h isotropic XY chain with random Lorentzian intersite coupling and transverse field at 
sites that depends linearly on the neighbouring couplings. We study in detail the wave 
vector- and frequency-dependent zz structure factor at different values of Hamiltonian 
parameters and temperature. We discuss the changes in the frequency profiles of the zz 
dynamic structure factor caused by the introduced correlated disorder. 
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Quantum spin chains with randomness have attracted a great deal of interest over the last 
decades. Apparently the simplest example of such a system is a random spin-1 XY chain the 
study of which is essentially simplified due to the fact that with the help of the Jordan- Wigner 
transformation it can be presented in terms of spinless noninteracting fermions. Starting from 
the early 70s several types of random spin-1 XY chains were discussed by using fermionization 
and Dyson's and Lloyd's models of disorder, as well as a numerical approach [1-7]. The interest 
in random spin-| XY chains has been strongly renewed in the 90s in view of the study of generic 
features of quantum phase transitions in disordered systems. As an example we mention here an 
exhaustive study of the transverse Ising chain by both the renormalization group and numerical 
means [8-10], and exact analytical and numerical treatment of XY chain [11-13]. 

In this paper, we deal with the spin-| isotropic XY chain with random Lorentzian exchange 
coupling and transverse field that depends linearly on the surrounding exchange couplings. Such 
a model has recently been examined in some detail [14]. The method developed by John and 
Schreiber [15] allowed one to calculate exactly the random-averaged density of magnon states 
for that model and, as a result, to study rigorously its thermodynamics. The most interesting 
property that the model with a correlated disorder exhibits is the appearance of a nonzero 
averaged transverse magnetization at the zero averaged transverse field [14,16,17]. This is 
conditioned by the change in the density of states due to the correlated disorder. Namely, the 
numbers of "magnons" with negative and positive energies at the zero averaged transverse field 
become not equal to each other. Unfortunately, the obtained analytical results pertain only to 
thermodynamics. The aim of the present paper is to study the effects of a correlated disorder on 
dynamics of transverse spin correlations, examining for this purpose the zz dynamic structure 
factor. To reveal the effects of the correlated off-diagonal and diagonal disorder it is necessary to 
analyse also the model with independent random Lorentzian exchange couplings and transverse 
fields. Such a study of dynamic properties requires the calculation of the zz time-dependent spin 
correlation functions and concerns the dynamics of the spin model conditioned by the exciting 
of only two magnons. Apparently, the evaluation even of the simplest zz time-dependent spin 
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correlation functions cannot be performed analytically, however, it can be done numerically, 
making use of the developed earlier finite-chain calculation scheme [18] (similar approaches are 
described in Refs. [9,10,19,20]). 

It should be stressed that models with the correlated disorder naturally arise while describing 
materials with the topological disorder. On the other hand, dynamic measurements are the 
basic experimental techniques in the study of such compounds. Although there are a few 
examples of real materials which are well described by a one-dimensional spin-| isotropic XY 
model (for instance, PrCl 3 [21]) the presented below exact numerical results do not pertain to 
any particular compound. Nevertheless, they still may be of much use for understanding the 
possible changes in observable quantities caused by the correlated disorder and can help to link 
the theoretical predictions and experimental data. 

Following Ref. 14 we consider a linear isotropic XY chain of N spins | in a transverse field 
governed by the Hamiltonian 
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where Q n is the transverse field at site n and J n is the exchange coupling between the sites n and 
n + 1. The J n are taken to be independent random variables with the Lorentzian probability 
distribution 

p(J n ) = —„ (2) 

1 ^ ( J n - J ) 2 + P K ' 

where Jo is the mean value of the exchange coupling and V is the width of distribution that 
controls the strength of the disorder. For the correlated off-diagonal and diagonal disorder the 
transverse fields are related to the intersite couplings and further it is assumed that 

n n -n = a I + I (3) 

where fi is the averaged transverse field at site [22]. One can easily find the probability 
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distribution for the random variable f2 



n 
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i.e. Vt n appears to be a Lorentzian random variable with the mean value fl and the width of 
the distribution \a\T. 

Mapping spin model (1) with the help of the Jordan- Wigner transformation onto spinless 

noninteracting fermions one obtains the Hamiltonian 

Y n N 
H — ^ ' firt ~\~ ^ ' c n A nm c m 

n=l n,m=l 



with 



that can be diagonalized by the transformation rj+ = J2f=i 9kjCj , = J2jLi9kjCj, where 

AT A^ N 

9kjAj m = Aj,gj, m , ~^^9kj9qj = $kqi ^ 9k j 9km = $jmi (5) 
j=l j=l k=l 

with the result if = £^ =1 A K (r?+?7 K - |) . Since = -±tp+tpj, where </?± = EjLi 9 P j (Vp ± , 
the calculation of the zz time-dependent spin correlation functions reduces to exploiting the 
Wick-Bloch-de Dominicis theorem with the outcome 

where the elementary contractions read 

AT / giApt e -iApt 

(V>j(t)ipt) = -{<Pj(t)<p~) = I 1 + e/3Ap + i^-flAp 

Af / e iApt e -iApt 

(<^ + (*Vm) = -(<Pj(t)V>t) = -J29pj9pm I 1 + e/3Ap - | ~ 



(7) 



Formulae (5) - (7) are the basic ones for the the presented below numerical study of the 
dynamic properties of spin model (1) - (4). For the given realization of random couplings (2) 
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and corresponding transverse fields (3) or for the given realizations of random couplings (2) 
and random transverse fields (4) one must first calculate the eigenvalues and eigenvectors of 
matrix A (5). Knowing its eigenvalues and eigenvectors one immediately obtains elementary 
contractions (7) and, therefore, the zz time-dependent spin correlation functions (6) that are 
directly related to the zz dynamic structure factor or susceptibility. Usually one is interested 
in the random-averaged quantities that come out as the result of averaging the computed 
quantities over many random realizations. The random-averaged quantities will be overlined. 
More details on the finite-chain calculation scheme can be found in Ref. [18]. 

In what follows we shall discuss the dynamics of transverse spin correlations in spin model 
(1) - (4), calculating for this purpose the transverse dynamic structure factor 
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The transverse dynamic structure factor (8) for a certain random realization with the help of 
Eqs. (6), (7) can be rewritten in the following form 

5(m + A p - Ag) 
n=o, ± i, ±2l ... P , g =i v . + e^)(l + e-/*,) 

[23]. For a uniform cyclic infinite chain (Q n = Q , J n = J ) the corresponding result reads 

6(u + A K > - A K '_ K ) 
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A K = Q + J cos k. 



(10) 



Evaluating the integral in Eq. (10) one gets the following expression for the transverse dynamic 
structure factor of the uniform cyclic infinite chain 

S K * |_ Js^k* i /sWk*-*"!! Fi ,a0\„*\fi ,o-/3A„*_ K N) ^ ^ < 2 | J Sin ^ | , 



S zz (k,lo) = < 



-JsinK*+Jsin(K*-K)| (l+e /3A K*)(l+e _/3A '«*- K ) ' 

0, otherwise 



(11) 
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Let us describe the sketched above numerical analysis of the dynamic properties of the 
considered spin system in more detail. To understand the accuracy of the numerical results we 
performed many additional calculations. In Figs, la, lb one can see the time dependence of 
the autocorrelation function {s Z N_(t)s Z N) — (s|) 2 fc> r N = 150 and N = 300, respectively. The 

2 2 2 

plots demonstrate the finite size effects that appear at t w 150 for N = 150 and t p=j 300 for 
N = 300. The same effects can be seen in the time dependence of (si5o(£) s 25o) ~~ ( s i5o) ( s 25o) 
computed for N = 300, which is depicted in Fig. lc. This correlation function appears with a 
delay at t p=* 100 and at times t p^ 200 exhibits a time behaviour influenced by the finite size of 
the chain considered. Figs. 2a and 2b demonstrate the difference between the sums 



(N = 300) with different n*. Because of increasing the time of delay in the appearance of 
correlation functions with large n, it is necessary to take a sufficiently large number of terms 
in the sum to reproduce correctly its time behaviour at large times (compare Figs. 2a and 
2b that correspond to n* = 50 and n* = 100, respectively). In practice one must restrict 
the calculations of the time- dependent spin correlation functions to some finite time of cut- 
off t c that generates wiggles in the frequency dependence of the dynamic structure factor. 
The wiggles can be removed by increasing the time of the cut-off to such values at which 
the spin correlations are already small enough or by increasing the value of e that smoothes 
the frequency profiles. Fig. 3 demonstrates how the numerical results approach the exact 
ones with increasing t c and decreasing e. Finally, let us add that acting in the described 
manner we reproduce the analytical results for the transverse dynamic susceptibility in a non- 
random case [24] (see also [25-27]). Figs. 4a - 4c show that in the random case one may take 
much smaller values of t c , since the random-averaged sum of correlation functions that yields 
S zz (k,,u>) decays essentially faster than in the non-random case (compare Figs. 4c and 2). 
Besides, with increasing the number of random realizations the resulting random-averaged sum 
of correlation functions becomes more regular (compare Figs. 4a, 4b and 4c). In our numerical 
calculations we considered spin-| transverse isotropic XY chains of N = 200, 300, 1800 spins 
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n=0,±l,...,±ra* 
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with J = — 1 (the results for S zz (n,u>) do not depend on the sign of exchange coupling), 
Q = 0,0.25,0.5,0.75,0.99 and T = 0, 0.1 at low {(3 = 1000) and high (f3 = 1) temperatures. 
We computed correlation functions (s Z N_(t)s Z N, ) — (s z n)(s z n, ) with n up to n* = 80, . . . , 800 for 

2 2 ' n 2 2 ' n 

the times up to t c = 150, . . . , 1200, put e = 0.005 and averaged the zz dynamic structure factor 
at least over 21000 random realizations. In the non-random case we used exact formulae (11), 
(12). The obtained within the frames of the described scheme results for the transverse dynamic 
structure factor of the non-random chain and the chains with a correlated and non-correlated 
Lorentzian disorder are presented in Figs. 5 and 6, respectively. 

Let us turn to the discussion of the obtained results. At first let us consider a non-random 
case (formulae (10) - (12), Fig. 5). To explain the observed frequency profiles one must take 
into account the fact that they reflect the dynamic properties of the magnetic chain conditioned 
by the exciting of two magnons with energies A K i = Q + J cos k' and A K « = Q + J cos n" 
for which oo = — A K > + A K » and k" — k' — k. Besides, at T = ((3 — > oo) the magnon energies 
must have the corresponding signs due to the Fermi factors involved into (10), namely A K > < 0, 
A K n > 0. Consider, for example, S zz (j,u) for Q = 0.5 (the curve for k — | in Fig. 5c and the 
dashed curve in Fig. 6b). Since A K = 0.5 — cos ft for T = 0, because of the Fermi factors, k' 
may vary only from — | to — ^ and, hence, k' — | varies from — 1| to — |, whereas for non-zero 
temperature k' varies from — n to n. The value of — A K > + A K /_| = —2 sin | sin (k' — |) is shown 
in Fig. 7a. Evidently, at T — 0, the lower frequency at which the non-zero value of S zz (^,uj) 
appears, Ui, is equal to cos ^ — cos | = 0.46592 . . . and the upper frequency uj U) after which 
S zz (j,cu) disappears, is equal to cos | — cos ^| = 0.75881 . . . ; for non-zero temperature u>i = 
and u u = 2 sin | = 0.76536 . . . (see Fig. 7a). The value of S zz (^, u>) is determined by the value 
of the slope of the curve depicted in Fig. 7a and this explains why S zz (^,ui) < S zz (j,uj u ) at 
T = 0, as well as why S zz (j,uj u ) — > oo at non-zero temperature. Note that, as it can be seen 
from Fig. 7a, the value of S zz (^,u) is determined by two pairs of magnons that satisfy the 
conditions uj = —A K i + A K » and k" = k' — k, however, at T = 0, because of the additional 
requirement A K / < 0, A K » > 0, only one pair of magnons can contribute to S zz (^,u). 
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Consider further, for example, S zz (^-, to) for f2 = 0.5 (the curve for k = in Fig. 5c and 
the dashed curve in Fig. 6d). For T — 0, k' may vary only from — | to |, whereas for non-zero 
temperature k' varies from — n to 7r. The value of — A K > + A k ,„2zl = —V% sin (k' — |) is shown 
in Fig. 7b. As it can be seen from this figure, for T = 0, as well as for non-zero temperature 
Ui — and uo u = = 1.73205 .... Besides, S zz (^f,uj u ) — > oo for any temperature. The abrupt 
change in S zz (^-,u) at T — at u> — 1.5 is due to the fact that for lower frequencies only one 
pair of magnons, because of the Fermi factors, contributes to S zz (^-,co), whereas for higher 
frequencies two pairs of magnons are involved in forming S zz (^-,u>). At non-zero temperature 
S zz , u>) is always conditioned by two pairs of magnons. 

Let us pass to random models (formula (9), Fig. 6). For such models the transverse dynamic 
structure factor is again conditioned by two magnons A p and A g , for which u> = — A p + A q and 
the quantity 

A(«,P,9)= X! elKn 9pj9 P ,j+n9qj9q,j+n 

n=0,±l,±2,... 

has a non-zero value. Besides, at T = 0, A p < 0, A q > 0. Apparently, there is no simple 
rigorous explanation for the behaviour of the random-averaged frequency profiles depicted in 
Fig. 6. However, it is easy to note that the kind of naive reasoning presented below does work 
for such a case. Consider at first the case k = j at low temperature. As it was shown above, 
the non-zero value of S zz (j,ui) in a non-random case was conditioned by two magnons with the 
energies 0.5 — cos ^ = —0.46592 . . . and 0.5 — cos | = 0, respectively. As it can be seen in Fig. 
8, where the random-averaged densities of magnon states p(E) = J2k=i ~ ^-k) obtained 
from (5) are depicted, such a pair of "magnons" does exsist for a = —1.01 (Fig. 8a) and does 
not exsist for a = 1.01 (Fig. 8b) or for the case of a non-correlated disorder (Fig. 8c). This 
observation is in agreement with the changes in the frequency profile due to different types of 
disorder shown in Fig. 6b (compare the curves 1, 4 and 2, 3 at u>i). In the non-random case the 
non-zero value of S zz (^, uj u ) was conditioned by two magnons with the energies 0.5 — cos | = 
and 0.5 — cos || = 0.75881 . . . , respectively. From Fig. 8 one can see that the density of 
states for such magnons is diminished because of the disorder that agrees with the changes in 
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S zz (j,uj u ) shown in Fig. 6b. 

Consider further the case k = ?f at low temperature. For the non-random case S zz (^f,uji) 
arises due to two magnons A| and A_z with the zero energy. As it can be seen in Fig. 
8, the disorder affects the density of states at E = for a = 1.01 and the non-correlated 
disorder more than for a = —1.01 that agrees with the changes in the frequency profile (Fig. 
6d). S zz (Jy,oj u ) is formed by two magnons with the energies 0.5 — cos | = —0.36602 . . . and 
0.5 — cos ^ = 1.36602 . . . , respectively. The density of magnon states for such energies is more 
diminished for the non-correlated disorder (Fig. 8c) than for the correlated one (Figs. 8a, 8b) 
which agrees with the smaller value of S zz (^-, uj u ) in the former case in comparison with the 
latter. 

To summarize, we extended the consideration of the spin-| isotropic XY chain with a 
correlated Lorentzian disorder presented in Ref. 14, examining numerically the dynamics of 
transverse spin correlations. We obtained the frequency dependences of the transverse dynamic 
structure factor at different values of the wave vector and temperature. We found the possible 
influences of the correlated disorder on the frequency profiles of the transverse dynamic structure 
factor. Within certain frequency regions the introducing of the correlated disorder may yield 
almost no changes in the value of S zz (k,cu), whereas the non-correlated disorder always erodes 
the frequency profiles of S zz (k,uj). The studied possible influences of the correlated disorder 
on the dynamic properties may be useful for the analysis of experimental data obtained in 
dynamic experiments for quasi-one- dimensional spin-| isotropic XY compounds. 

The authors are grateful to Prof. M. Shovgenyuk for providing the possibility to perform 
the numerical calculations. The paper was discussed at Magdeburg University and Dortmund 
University. O. D. is grateful to Prof. J. Richter and Prof. J. Stolze for their warm hospitality. 
He is also indebted to Mrs. Olga Syska for the financial support. 
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List of figure captions 
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FIG. I. Time dependence of the transverse spin correlation functions r J J+n (t) = Re {sj(t)sj 
-( s j)( s j+n)] for J o = -1, = 0.5, r = at (3 = 1000. a) N = 150, j = j + n = 75; b) 
N = 300, j = j + n = 150; c) N = 300, j = 150, j + n = 250. In this figure and in Figs. 
2, 4 we plotted only the real parts of correlation functions since their imaginary parts exhibit 
qualitatively the same behaviour. 



FIG. 2. Time dependence of K?(t) = ReE n=0 ,±i,...,±„.e lra [{s z 150 (t)s z 50+n ) - (sf 50 ) (s z 150+n ) 
k = tt for Jq = -1, tt = 0.5, r = 0, iV = 300 at (5 = 1000 for n* = 50 (a) and n* = 100 (b). 



FIG. 3. S zz (k,u) at re = f (1), k=\ (2), « = f (3), « = ^ (4), « = ?r (5) for the uniform 
chain with J = —1, Q = 0.5 at j3 — 1000: exact results (11), (12) (dashed curves) versus 
numerical ones (solid curves), a) N = 280, j = 140, rf = 100, t c = 200, e = 0.01; b) N = 1800, 
j = 900, n* = 800, i c = 1200, e = 0.01; c) = 1800, j = 900, n* = 800, t c = 1200, e = 0.005. 



FIG. 4. Time dependence of Rf(t), k = n for J = -1, fi = 0.5, T = 0.1, a = 1.01, N = 300, 
n* = 100 at (3 — 1000; the random-averaged quantity comes as the result of averaging over 
1000 realizations (a), 11000 realizations (b), and 21000 realizations (c). 



FIG. 5. Frequency dependence of the transverse dynamic structure factor (11), (12) for J = —1 
and different values of transverse field f2 = (a), f2 = 0.25 (b), f2 = 0.5 (c), f2 = 0.75 (d), 
f2 = 0.99 (e) at different values of wave vector k = | , | , | , ^, tt (from left to right) at 
P = 1000. At high temperatures the frequency profiles of S zz (k, uj) are the same for all values 
offio (f). 
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FIG. 6. Frequency dependence of the random-averaged transverse dynamic structure factor (8) 
at different values of wave vector k — | (a), re = | (b), k — f (c), k — (d), k — (e), 
k = 7T (f) for model (1) with J = -1, Vl = 0.5, r = 0.1 at f3 = 1000. 1) correlated disorder (3) 
with a = —1.01; 2) correlated disorder (3) with a = 1.01; 3) independent exchange couplings 
and transverse fields, the latter are distributed according to probability distribution (4) with 
|a|T = 0.101; 4) non-random case T = (dashed curves). 

FIG. 7. — A K i + A K /_ K versus k' for the non-random chain with J = —1, f2 = 0.5. a) k — ^; 
b)« = f. 



FIG. 8. Density of states p(E) for model (1) with J = — 1, fio = 0.5, T = 0.1. a) correlated 
disorder (3) with T = 0.1, a = —1.01; b) correlated disorder (3) with T = 0.1, a = 1.01; c) 
independent exchange couplings (2) with T = 0.1 and transverse fields (4) with \a\T = 0.101; 
the density of states for the non-random case T = is depicted by dashed curves. 
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Figure 1: Time dependence of the transverse spin correlation functions rjj +n (t) 



Re [<sj(t)sj +n ) -(s*)(s* +n )J for J = -1, fi = 0.5, r = at /? = 1000. a) iV = 150, 
j = j + n = 75; b) N = 300, j = j + n = 150; c) N = 300, j = 150, j + n = 250. In this figure 
and in Figs. 2, 4 we plotted only the real parts of correlation functions since their imaginary 
parts exhibit qualitatively the same behaviour. 
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Figure 2: Time dependence of Rf (t) = Re E„=o,±i,...,±„* e lKn [{s z 150 (t)s z 150+n ) - (s* 50 ) (s z 150+n ) 
k = 7T for J = -1, tt = 0.5, r = 0, N = 300 at (3 = 1000 for n* = 50 (a) and n* = 100 (b). 
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Figure 3: S zz (k,uj) at re = f (1), re = f (2), re = ^ (3), re = f (4), re = vr (5) for the uniform 
chain with J = — 1, f2o = 0.5 at (3 — 1000: exact results (11), (12) (dashed curves) versus 
numerical ones (solid curves), a) iV = 280, j = 140, n* = 100, t c = 200, e = 0.01; b) N = 1800, 
j = 900, n* = 800, t c = 1200, e = 0.01; c) = 1800, j = 900, n* = 800, t c = 1200, e = 0.005. 
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Figure 4: Time dependence of B%*(t), k = tt for J = -1, fi = 0.5, T = 0.1, a = 1.01, 
N = 300, n* = 100 at (3 — 1000; the random-averaged quantity comes as the result of averaging 
over 1000 realizations (a), 11000 realizations (b), and 21000 realizations (c). 
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Figure 5: Frequency dependence of the transverse dynamic structure factor (11), (12) for 
J = —1 and different values of transverse field f2 = (a), f2 = 0.25 (b), f2 = 0.5 (c), 
tt = 0.75 (d), Q = 0.99 (e) at different values 8 of wave vector « = f , f , § , ^, 2p, tt (from left 
to right) at /? = 1000. At high temperatures the frequency profiles of S zz (k,u>) are the same 
for all values of fi (f)- 




Figure 6: Frequency dependence of the random-averaged transverse dynamic structure factor 

(8) at different values of wave vector n — | (a), k — | (b), k — f (c), k — (d), k = (e), 

« = 7r (f) for model (1) with J = -1, Vt = 0.5, r = 0.1 at f3 = 1000. 1) correlated disorder (3) 
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with a = —1.01; 2) correlated disorder (3) with a = 1.01; 3) independent exchange couplings 
and transverse fields, the latter are distributed according to probability distribution (4) with 
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Figure 8: Density of states p(E) for model (1) with J = — 1, f2 = 0.5, T = 0.1. a) correlated 

disorder (3) with T = 0.1, a = —1.01; b) correlated disorder (3) with T = 0.1, a = 1.01; c) 

independent exchange couplings (2) with T = 0.1 and transverse fields (4) with \a\T = 0.101; 

the density of states for the non-random case T = is depicted by dashed curves. 
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